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Summary 

Work  completed  under  this  grant  on  polar  fluids  includes  an  extension 
of  perturbation  theory  to  polar  molecules  with  point  octopoles  and 
hexadecapol es  as  well  as  a  preliminary  investigation  of  molecules  with 
non-spherical  cores.  In  addition,  the  ordering  of  dipolar  molecules  near 
a  flat  electrified  wall  has  been  investigated  in  the  linearized  hypernetted 
chain  (LHNC)  and  quadratic  hypernetted  chain  (QHNC)  approximations.  This 
has  led  to  a  new  study  of  el ectrostriction  and  non-linear  effects  in  the 
polarization  density  of  dipolar  fluids. 


This  report  describes  research  on  polar  fluids  carried  out  at  the 
University  of  Maine  at  Orono  under  contract  from  the  Office  of  Naval  Research 
Apart  from  the  principal  investigator  (Professor  J.  C.  Rasaiah),  a  visiting 
fellow.  Professor  Noel  Thompson  from  Australia,  a  graduate  student,  Mr.  John 
Eggebrecht  and  a  postdoctoral  fellow.  Dr.  Dennis  Isbister  from  the  Royal 
Military  College,  University  of  New  South  Wales  in  Australia,  worked  on  this 
project;  the  last  two  receiving  partial  financial  support  under  the  present 
contract. 

The  primary  objective  of  this  investigation  was  to  extend  the  theory 
of  polar  fluids  developed  by  Stell,  Rasaiah  and  Narang^  to  (a)  spherical 
molecules  with  point-octopoles  and  point-hexadecapol es  which  could  serve 
as  a  reasonable  model  for  methane  and  sulfur  hexafluoride;  (b)  non-spherical 
molecules  with  dipoles  and/or  quadrupoles  such  as  HCN  and  O^.  A  second 
objective  was  to  extend  our  study  of  dipoles  to  surface  adsorption  near 
an  electrified  wall . 

A  Monte  Carlo  calculation  of  the  free  energy  changes  that  accompany 
the  switching  on  of  polar  interactions  between  the  molecules  was  planned 
as  an  unambiguous  test  of  the  accuracy  of  the  theory  when  applied  to  non- 
spherical  cores.  Some  aspects  of  this  work  are  still  in  progress,  while 
parts  (a)  and  (b)  of  the  research  described  have  been  completed.  In 
addition,  a  very  promising  and  important  theoretical  investigation  of  the 
adsorption  of  dipoles  at  a  wall  in  the  presence  of  an  electric  field  has 
been  initiated  with  the  assistance  of  Dr.  Dennis  Isbister  and  John 
Eggebrecht.  In  what  follows,  the  research  that  has  been  completed  with 
support  from  ONR  is  described  in  greater  detail  and  future  plans  to  continue 
this  work  are  outlined. 
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I.  Perturbation  Theory  for  Polar  Fluids 

While  the  potential  energies  of  interaction  amongst  point-dipoles 

and  point-guadrupol es are  wel 1  known,  there  are  inconsistencies  and  ambigu- 
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ities  in  the  literature  when  higher  order  multipoles  were  included. 

It  was  therefore  decided  to  calculate  these  energies  from  first  principles, 
limiting  our  attention  only  to  octopolar  and  hexadecapolar  interactions 
beyond  quadrupolar  charge  distributions.  It  became  necessary  therefore  to 
perform  the  tedious  contractions  of  multipolar  and  field  tensors  the 
eight  index  field  tensor  for  the  hexadecapole-hexadecapole  energy  contains 
nearly  700  terms.^  The  two-body  potentials  beyond  the  simple  dipole-dipole, 
dipole-quadrupole  and  quadrupole-quadrupole  interactions  were  found  to  be 
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where  i^,  e,  u,  and  t  are  the  dipole,  quadrupole,  octopole  and  hexadecapole 
moments  and 

COS0..  =  cose. cose  .-sine. sine. cos(()). -41.) 

I J  '  J  *  J  ’  J 

where  (9^,4)^)  are  the  polar  and  azimuthal  angles  of  orientation  of  the 
axis  of  symmetry  of  the  ith  molecule. 

The  potential  energy  functions  were  then  expanded  in  spherical 
harmonics,  to  allow  analytic  evaluation  of  the  angular  integrals  in  the 
x-expansion  of  the  free  energy 
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where  and  ((x.^)  are  the  spherical  harmonics,  the 

spatial  and  angular  variables  of  the  potential  energy  can  be  separated. 

Using  the  orthogonality  properties  of  the  spherical  harmonics,  and  the 
previously  derived  potential  energy  functions  we  have  found  the  following 
expansion  coefficients  for  the  higher  order  multipole  interactions  (t-j  =  1, 
2,  3,  and  4  correspond  to  dipole,  quadrupole,  octopole  and  hexadecapole 
moments) 
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The  general  pattern  of  these  coefficients  allows  a  check  on  the  tedious 
and  mistake-prone  calculations. 
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With  these  results  the  terms  of  0(a  )  in  the  free  energy  expansion  can 
be  calculated.  To  go  beyond  this,  using  a  Pade  approximant  for  the  free 
energy  in  the  manner  proposed  by  Stell,  Rasaiah  and  Narang,  requires  at 
least  the  term  of  0(x  ).  To  this  end  the  corresponding  angularly  averaged 
three-body  potentials  for  higher  multipoles  were  needed.  They  were  not 
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available  in  the  literature  but  were  derived  in  the  manner  of  R.  J.  Bell 
and  Rasaiah  and  Stell.^ 


5 


Using  quantum-mechanical  perturbation  theory  Bell  has  shown  that  the 
non- additive  three-body  potential  due  to  fluctuating  dipoles  and  multipoles 
may  be  expressed  as  a  function  of  the  intermolecular  distances  and  the 
interior  angles  of  the  triangle  which  they  form 
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Rasaiah  and  Stell  recognized  that  the  corresponding  three-body  term  in  the 
perturbation  theory  of  polar  fluids  had  the  same  form 
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The  geometrio  terms  W  ,  W  ,  W  ,  and  W  have  yet  to  be  obtained. 
Work  on  this  is  prooeeding, 
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To  determine  the  speoifio  ooeffioient  Z  1  2  3  in  the  three-body 

potential  we  made  use  of  the  spherical  harmonic  expansion  coefficients  and 
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the  method  outlined  by  Rasaiah  and  Stell  to  obtain 
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From  these  sequences  one  may  predict  the  values  of  Z  1  23  even  when  the 
geometric  terms  beyond  W224  have  not  yet  been  obtained.  We  expect  that 
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To  complete  this  portion  of  the  work  we  needed  to  evaluate  the  remaining 
four  geometric  factors  to  obtain  the  effective  angularly-averaged  three-body 
potentials  which  contribute  to  the  term  of  0{ <  )  for  all  binary  mixtures 
with  multi  pole  moments  through  the  hexadecapole  moment.  The  integrals  of 
the  form 


123^'"l2’'"l3’''23 
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were  then  numerically  evaluated  using  a  method  previously  developed  by 
Rasaiah,  Larsen  and  Stell.^  The  result  is  the  first  and  second  terms  for 
the  lambda  expansion  from  which  a  Pade  approximant  for  the  free  energy 
can  be  formed. 

In  order  to  extend  these  perturbation  methods  to  systems  with  non- 
spherical  cores  we  must  first  properly  characterize  the  reference  system. 

To  that  end,  we  proceeded  with  a  computer  simulation  of  linear  non- 
spherical  molecules.  The  objective  was  to  compute  an  ensemble  averaged 
angularly  dependent  distribution  function  in  terms  of  the  coefficients  in  a 
spherical  harmonic  expansion  of  the  distribution  function 
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From  these  coefficients  one  may  obtain,  by  the  orthogonality  of  the  spher¬ 
ical  harmonics  in  the  expansions  of  the  distribution  and  potential  func¬ 
tions,  a  very  simple  expression  for  the  term  of  0('‘)  for  a  linear  molecule 
with  a  quadrupol e-moment 
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where  is  the  diameter  of  one  of  the  atoms  in  the  diatomic  'i.Dlec^ie  and  . 
is  the  reduced  density. 

Monte  Carlo  simulations  of  the  equilibrium  properties  of  ^  nomonuclear 
diatoicic  molecule  with  molecular  dimensions  characteristic  of  tl-  were 
pei'Tc ',iec  to  generate  an  accurate  set  of  coefficients  g  [r ,  apoearing 

j  i  '  ')]• 

I 

ir.  tne  spnerical  harmonic  expansion  oi  tr,e  radial  di str i rtu i on  'unctun 
0(1',  ,  Nearly  a  million  configurations  were  generated  ir  ca  1  c  a  1 1  nq 

thf  results  deoicted  in  Fig.  1.  T'-.e  moi'r-nts. 


wnicn  appear  in  the  term  of  0('',)  in  the  rree  energy  decreased  steadMy  by 
cnlyi'.S'  over  the  last  600K  configurations.  These  calculations  are 

protariv  "ore  accurate  than  earlier  detenni nations^ of  g  (r).  even 

'  1  ■  C'’’ 

tno.^cri  tne  programs  used  were  essentiallv  the  same.  Tnis  is  primarily 
betaiise  a  larger  number  of  configurations,  with  more  f  recer  t  sar.:  linos  at 
s::..  Ic  intervals,  were  generated.  Ir  Ta.''’e  I,  the  newer  estimates  tne 
ton  e*  ■)  are  compared  with  older  data,  for  reduced  auadrupo'.e  '  orents 
coriespona i ng  the  critical  point  and  melt'ing  point  of  CI2  at  a  reduced 
density  of  0.5. 

Table  I 

Comparison  of  the  term  of  0(x)  obtained  by  numerical  integration 
of  (7)  using  the  data  for  g^oglr),  q^.^,(r),  and  from  the 

independent  studies  of  Street  and  Tifdesley  (S&T;  afid  Eggebrecht 
and  Rasaiah  (E&R) 
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.283 
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.493 

.571 
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.723 

i  .975 

.686 

.798 

.907 

1.006 
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A  number  of  difficulties  became  apparent  from  tbe  discrepancies  in  the  table. 
The  sources  of  error  in  this  calculation  were  found  to  be:  1)  the  large 
fluctuations  in  the  moments  of  sensitivity  of  the  function 

to  the  method  of  integration.  The  greater  similarities  in  our  estimates  of 
the  integral  may  be  attributed  to  the  smaller  increments  used  in  our  version 
for  the  radial  distance. 

II.  Adsorption  of  Dipoles  and  Ions  at  a  Wall  in  the  Presence  of  an  Electric 
Fi  eld 

While  our  understanding  of  the  behavior  of  liquids  has  improved 
dramatically  over  the  past  ten  years,  the  physics  and  chemistry  of  surface 
phenomena  is  still  at  a  primitive  stage.  One  of  the  outstanding  problems 
in  this  area  is  the  elucidation  of  the  density  profiles  of  ions  and  dipoles 
at  a  charged  surface;  a  detailed  knowledge  of  this  is  fundamental  to  our 
understanding  of  phenomena  ranging  from  the  behavior  of  electrodes  in 
batteries  to  selective  adsorption  and  transport  across  membranes.  It  is 
convenient  to  study  ions  and  dipoles  at  a  surface  separately  before  trying 
to  gain  insight  into  the  behavior  of  a  mixture  of  these  two.  We  have, 
therefore,  embarked  upon  a  program  in  which  the  following  problems  will  be 
studies  in  sequence. 

1)  Dipoles  at  a  wall  in  the  presence  of  an  electric  field. 

2)  Ions  at  a  charged  interface. 

3)  Dipoles  and  charges  at  a  wall  in  the  presence  of  an  electric  field. 

Investigation  of  the  first  problem  has  progressed  very  rapidly  in 

collaboration  with  Dr.  Dennis  Isbister  and  John  Eggebrecht.  Professor  Noel 
Thompson,  a  visitor  on  sabbatical  leave  from  Australia,  Dr.  Dennis  Isbister, 
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and  I  have  made  progress  on  the  second  problem,  while  work  on  mixtures  of 
ions  and  dipoles  at  a  wall  has  only  just  begun.  In  what  follows,  we  provide 
a  brief  description  of  our  work  on  dipoles  at  a  wall. 

The  effect  of  a  wall  is  obtained  by  allowing  a  mixed  system  of  at  least 
two  species  1  and  2  to  change  in  such  a  way  that  the  density  of  one  species 
(2  say)  goes  to  zero  accompanied  by  an  increased  in  diameter  R^.  It  is  the 
growth  of  this  single  particle  (2)  that  emerges  finally  as  a  wall,  and  if 
the  wal 1 -particl e  had  a  dipole-moment,  the  zero  density  -  infinite  radius 
limit  could  be  made  to  produce  an  electric  field  as  well  emerging  at  any 
chosen  angle  from  the  wall.  The  density  profile  of  each  fluid  species 
near  the  wall  is  given  by 

(r)  =  lim  lim  p  ''[h^-,  (p  ,;l!2)+l  ] 

where  o^^is  the  bulk  density  of  species  1,  and  ^21^'^^  total  wal  1 -parti cl e 
correlation  function  which  has  the  usual  invariant  expansion. 

Our  initial  investigation  of  dipoles  at  a  wall  used  the  linearized 
hypernetted  chain  (LHNC)  closure  for  the  wal 1-particle  and  particle-particle 
interactions . 

It  is  known  that  this  approximation  leads  to  fairly  accurately  bulk 
properties  for  dipoles.  We  found  that  the  contact  values  for  the  wall- 
particle  correlation  function  are  much  higher  than  those  obtained  in  the 
mean  spherical  (MS)  approximation.^ 

The  linearized  hypernetted  chain  closure  is  obtained  by  taking  the 
first  term  in  the  expansion  of  the  angular  part  of  ln(h2i+l)  which  appears 
in  the  hypernetted  chain  (HNC)  approximation  for  the  direct  correlation 
function  C^i .  If  the  second  term  is  retained  as  well,  a  more  accurate 
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description  of  the  system  results.  This  approximation,  known  as  the 

8  9 

quadratic  hypernetted  chain,  QHNC,  ’  reveals  a  new  phenomenon  that  is 
absent  in  the  MS  and  LHNC  approximations.  This  is  the  change  in  the 
density  of  the  fluid  in  an  open  system  when  the  electric  field  is  turned 
on,  and  is  called  el ectrostriction .  The  leading  term  in  the  relative 
change  in  density  at  an  infinite  distance  from  the  wall  was  found,  in 
the  QHNC  approximation,  to  be  given  by  a  term  of  0(E  ) 


(2)  ^  Ap_ 
h  p^o 


(e-1) 


2  E 


2 


Q 


where  E  is  the  electric  field,  Q  the  inverse  compressibility,  l  the  di- 

4-  2  0 

electric  constant,  and  y  =  y-  p.|  6,  with  p  =  1/kT.  Approximations 
beyond  the  QHNC  theory  systematically  generate  term  of  higher  order,  so 
that  this  approximation  must  contain  the  complete  electrostriction  term 
of  0(E  )  in  the  HNC  approximation.  However,  a  discrepancy  in  was 
found  on  comparing  our  result  with  the  thermodynamics  of  electrostriction 
in  an  open  system. Since  the  HNC  approximation  and  all  of  the  theories 
derived  from  it  ignore  bridge  diagrams,  the  discrepancy  in  the  term  of 
0(E  )  between  the  QHNC  result  and  thermodynamics  must  come  from  neglect  of 
these  diagrams.  The  lowest  order  bridge  diagram  was  then  calculated  and 
found  to  contribute  significantly  to  electrostriction.  Adding  this  to  the 
QHNC  result  of  0(E  ),  we  obtained  (resolving  the  above  discrepancy) 


K  _r  B(c-l)^  5(.-l)^6y1  E^ 

h  ^24Ttp.|Uy  128:tp^o  Jq 

which  shows  that  the  bridge  diagrams  cannot  be  ignored  in  calculations  of 
the  local  density  at  an  infinite  distance  away  from  the  wall.  It  appears 
likely  that  they  affect  the  density  profiles  also  when  the  dipoles  are 


13 


closer  to  the  wall.  This  is  an  important  result,  because  all  of  the 
theoretical  studies  of  the  ordering  of  dipoles  and  ions  at  a  wall  that 
have  been  carried  out  so  far  ignore  these  bridge  diagrams. 

Our  study  of  the  LHNC  approximation  for  the  wall -parti cl e  distri¬ 
bution  function  leads  to  the  constitutive  relation 

^{°°)  ~  (e-1)  ?{“>)/4Tr 

between  the  polarization  denisty  ?{“>)  and  the  electric  field  f  with  the 
same  dielectric  constant  c  as  that  obtained  in  relating  c  to  the  two- 
particle  correlation  function  of  a  dipolar  fluid  in  the  absence  of  an 
electric  field.  The  QHNC  approximation,  however,  yields  in  addition 
non-linear  terms  in  the  electric  field  for  the  polarization  density. 

We  plan  to  extend  this  work  to  mixtures  of  dipoles  in  order  to  study  the 
selective  adsorption  of  dipoles  near  an  electrified  wall. 
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Fig.l  The  coefficients  of  the  spherical  harmonic  expansion  of  the  angle  dependant 
distribution  function  after  two  million  configurations. 
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A  theory  of  elect  rostriction  which  tollow'*  from  studies  of  dipolar  I'lderinj:  .ii  an  electrified  wall  is  discussed  in  the  qua¬ 
dratic  hypernetted  chain  approviniaiion.  Bridge  diagrams  for  the  wall  particle  correlation  functions  contribute  significant¬ 
ly  to  electrostriction  even  to  lowest  order  in  the  electric  field.  The  form  of  the  constitutive  relation  between  the  polariza¬ 
tion  density  and  the  field  in  strong  tields  is  discussed. 


I .  Introduction 

The  ordering  oT  dipoles  at  a  wall  Irotti  w  hich  an 
electric  field  emerges  has  been  si  tidied  in  the  mean 
spherical  (MS)  [1 1  and  linearized  hvpernettcd  chain 
(LHNC)  12]  appro.ximations.  Wc  leport  here  that  elec¬ 
trostriction  appears  as  an  added  lealuic  when  the  the¬ 
ory  is  carried  beyond  the  LHNC  appioMinatum  for 
the  wall  particle  correlation  functions,  leading  to  a 
molecular  theory  of  this  phenoinenon  in  w  hich  graphi¬ 
cal  analysis  and  integral  equation  approMinai ions  that 
are  ubiquitous  in  the  theory  of  lluids  may  be  exploit 
ed.  Of  the  hierarchy  of  approximations  geneiatcd  by 
the  hypernetted  chain  equation  (lINCi.  we  find  that 
the  quadratic  hypernetted  chain  (01  INC)  approxima¬ 
tion  is  the  first  to  predict  electrostriction.  The  leading 
term  in  the  relative  change  in  density  an  infiitite  dis¬ 
tance  away  from  a  flat  wall  is  of  0(/  *  I,  where  H  is  the 
magnitude  of  the  local  field.  Approximations  beyond 
the  QHNC  theory  systematically  generate  terms  of 
higher  order,  so  that  the  QHNC  approximation  also 
contains  within  it  the  complete  electrostriction  term 
of  0(A'^)  in  the  HNC  approximation. 


.Ml  of  '.he  ilicories  ot'  electrostriction  which  have 
their  genevtv  in  the  IlNf  .  i  proximation  ignore  the 
bridge  diagrams,  hut  we  find,  by  comparing  the  QIINC 
theoiy  with  the  thermodynamics  of  electrostriction 
|.’!.  that  these  diagrams  must  contribute  significantly 
even  to  lowest  order  in  /  This  leads  us  to  believe  that 
they  play  an  equally  important  role  in  determining  the 
local  density  when  the  dipoles  are  closer  to  the  wall. 

lepoi!  hete  the  haiest  outline  of  our  molecular 
theoiv  lot  an  open  system  |4|  in  which  contributions 
liotii  the  OlINt  approximation  to  0(/-.'^)  and  from 
the  bridge  diagrams  to  the  same  order  in  the  electric 
lield  and  to  lowest  otdei  in  the  density  and  the  dipole 
moineni  aie  evaluated  analytically.  Our  theory  can  al¬ 
so  be  extended  to  highei  order  in  /;  but  this  will  not 
be  pursued  heie.  The  thermodynamic  discussion  of 
Kirkwood  and  Oppenhcim  |.T|  however  is  limited  to 
deriving  the  elect  tost  net  ion  term  of  0(/f*).  since  it  is 
based  on  a  thermodynamic  perturbation  theory  of 
first  order  in  /.  ^.  that  uses  as  reference-state  input  the 
linear  constitutive  relation  between  the  polarization 
density  />(“")  and  the  local  electric  field  E  for  a  field- 
independent  e. 
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/'(»=)  =  (e  1 )  £’(oo)/47r.  ill 

I  111  I  hoi  (Jeiails  ufour  work  and  exierisions  i<>  n 
will  he  pieseiilcd  elsewliero  (4|  but  we  remark  tiuii 
oiii  stiid\  yields  the  eonstitutive  relation  with  the 
same  value  ot  e  as  that  obtained  from  from  lelaline  < 
to  the  two-particle  orientational  correlation  function 
in  the  absence  of  a  field,  and  also  provides  a  basis  tot 
deiemiimiij;  non-linear  terms  in  the  electric  Held  foi 
the  polarization  density. 

2.  Eleetrostriction  in  an  open  system 

We  employ  the  technique  discussed  in  refs.  1 1 .2|  to 
iienerate  an  electric  field  by  taking  a  mixture  ol  dipo¬ 
lar  hard  spheres  to  the  zero-density  limit  of  species  2 
with  an  attendant  increase  in  the  radius  i  to  infinii> . 
under  the  constraint  that  its  dipole  moment  wi  divid¬ 
ed  by  the  cube  of  the  radius  of  the  excluded  volume 
A’  ■),  =  A  ■)  +  A I  is  a  constant  A.'f, .  The  elect  ric  field  /  , 
emerging  from  the  wall  is  related  to  b>  [1.2) 

A  1  =  cos-d2  11*  *■* 

wheier'i  is  a  unit  vector  dependent  upon  On  ll  .2|. 
the  angle  which  the  wall  dipole  makes  with  the  wall 
luirmal  //  Note  that  A  i  is  constant  when  the  wall  di 
pole  orienlation  Cj  is  fixed.  The  relationship  beiweci' 

A  1  and  the  Maxw’ell  field  A' m  the  fiuid  is  of  the  tonn 

A  =  |.’ t2e  +  1 1)  AN  +  terms  non-hnear  in  A.  T.  t  'l 

lleie  the  non-linear  terms  aie  ol  inagnilude  const  /  2 
+  .  .when  di  ~  0.  ^nd  e  is  the  dielectric  constani  ol 
the  Huid  at  zero  field.  In  the  notation  of  ref.  |2| .  the 
wall  particle  eorrclation  function  has  the  nivaiiam 
expansion 

/i-ijlr.A.T.Oi )  “  //2|fc)  +  /i2ilc)/fl2.  1 ) 

+  /;^l(rlA(:,l) +.  ..  (4| 

where  Al2, 1 )  =  S2  •*]  and  /A(2, 1 1  =  S2’(-^nn  U  I  Si 
in  which  Si  and4|  ate  unit  vectors  in  the  directions  ot 
the  wall  dipole  and  fluid  dipole  tespeetivelv .  r  is  the 
distance  ofa  fluid  dipole  (of  orientation  fl|  I  ftom  the 
wall,  and  U  is  the  unit  tensor.  In  this  wall  limit.  AAi 

-  *  OO^ 


where  is  short -langed  and  AN;  is  a  constant  re¬ 

lated  to  the  electiic  field  AN  isce  'Cclion  2  ).  Hence 

lim  /iij (cl  =  .' AN|  IM 

file  oiigiii  I'l  elect  i‘  'vi  I  let  loll  in  our  theo¬ 
ry  is  the  coupling  that  can  exi-.i  between  liQ^  (c)  and 
/iij  (r);  this  feedback  begins  w  itli  the  01  IN(  appiox- 
imation 

We  define  the  electrosliiction  elte,i  A;,  as  the  rela¬ 
tive  change  in  density  ot  ilie  bulk  H.iul  when  the  elec¬ 
tric  field  IS  sw  itched  on.  so  iliai 

Ay,  ^  Ap/p*,' =  AN  1,  ("1 

where 

)dn,  (h) 

=  477  for  dipolai  hard  spbetes  and  p'|'  is  the  density 
of  the  bulk  fluid  when  A  s  -  0  It  A  is  the  correspond¬ 
ing  asymptoiie  liniil  ot'ihe  angiilaily  .iveiaged  diiecl 
corielation  tunclion  (’21 1 ,11^  l.we  tind  |4|  that. 

Ill  the  absence  ol  inoleculai  polai  i/.d'ihi \ 

A/,  =A.  (A.  <'•» 

where  (A  is  the  inveise  compiessihihiy  >'•  the  lluid. 

0  =■■  I  (p*,’'U-i  /  ciiif.n,  .n  N  df2,  dn,  sir  lUti 

and  c'jilf.  n,  .  ft-i  I  IS  the  diie^i  ^.nte’.ilion  function 
of  the  bulk  fluid.  The  ihennodv  n.iini,  theory  ol  elec- 
trostiiction  |.'  1  loi  an  open  s\  sieni  gives  A ,,  toOiA  ’  I  as 

A/*’*  =  (|3-'S7t)((1c  dp'i'i/  -  Ill) 

where  P~  1  /A  T.  A  is  Itolt.-niann's  .onsl.nit .  7  is  the 
absolute  temperal urc  and  the  suj'ei s, ; ipi  ( 2 )  means 
the  teim  of  0(/  * ). 

-T  The  OHNf  approximation  for  eleetrostriction 

By  considering  the  asy  niptotw  Inn  it  ol  i^jlc)  m  the 
invariant  expansion  li'i  the  wall  ivailieie  direct  corre¬ 
lation  (unction  ('21  (-■-  ^  ;  ^^1  •  "-e  fi'id  (4|  that  ANj  is 
telated  to  A  ->  in  the  0IIN(  appioximal ion  by 


O  t  -  \  -  I- 
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;3»i,/;;;(3cos202  =  AS, |2(;?,,(2A:i,p*/«fl ) 

+  (9  (  A^iiP?«ll)  3A,,/(1  +A,,)1,  (12) 

whole  the  Q.(  )  functions  are  defined  in  ref.  [2 1 , 

and  the  constant  K  n  is  related  quite  itenerally  (i.e.  not 
lioi  111  the  QMNC  theory)  to  the  lluid  dipole  moment 
III,  by  [4  6 1 

3v  =  c>.(2Ai,p‘i'/?|i)  A,ip‘,’a],).  (13) 

wheie  V  =  Anm\p^0j^) .  The  dielectric  constant  e  of 
the  fluid  IS  also  quite  generally  given  hy  |(i| 

e  =  (.M2AnP?Ajii  Ai,p‘;A-Ji).  (14i 

We  also  find  from  the  r  -♦  <»  limit  of  the  OHN(  ap- 
pioxiniation  for  c'2i 

a;.-' =1  A5,(3  cos-Oj  +  I)  (15) 

and  on  using  (2),  (3),  (<)')  and  ( 12)  -  ( 14)  we  have 

a;,-  ’  =  ( /3' 24jrp°  v)  (e  -  1 )-  A"  /Q.  ( 1  h  I 

where  the  superscript  (2)  means  again  the  term  of 
C’>(/  -  ).  Note  that  in  (16)  A/,"’  is  independent  ol  the 
inclination  6-,  of  the  wall  dipole,  Consistency  with  the 
thermodynamic  theory  (11)  demands  that 

p',' i)e  9p',' =(e  l)r'3.f.  (IT) 

This  relation  is  indeed  satisfied  by  the  simple  Debye 
etiiijiion  (e  -  1  )/(e  +  2)  =,v.  but  the  more  exact  ex¬ 
pression  [7,81  (for  arbitrary  y  and  p*^  -*0) 

ic  l)(e  +  2)  =  y  H.i-W  ..  (18) 

introduces  a  discrepancy  of  0(>'^)in(17).  Since  higlier 
order  approximations  derived  from  the  I1N(  equation 
do  not  produce  additional  contributions  to  A;,  ol 
0(7  -  ),  the  apparent  inconsistency  lies  in  ignoring  the 
bridge  diagrams. 


4.  Bridge  diagrams  of  O(A^) 


The  bridge  diagram  of  OtA^)  and  to  lowest  order 
in  the  fluid  density  can  be  represented  graphically, 
when  :  by 


A*(2.1)  = 


(l‘») 


where  < - -  -  =  TAii  Z)(2.  r  ).  and 

-)  ^ 


is  the  three-particle  direct  correlation  function  fylAj . 
.V3.  A'^)  for  the  fluid  with  A',  =  (r,,  |,  A*(2.  1 )  con¬ 

tributes  to  the  wall  -  particle  direct  correlation  func¬ 
tion.  In  (19)  the  open  circle  2  is  the  wall  dipole  root 
point,  the  half  black  circle  »  signifies  angular  integra¬ 
tion  over  the  orientations  of  the  bulk  fluid  dipole  I 
and  the  black  circles  (field  points)  lepresent  spatial 
and  angular  integrations  id  fluid  dipoles.  To  evaluate 
(19)  to  lowest  order  in  p^j’  and  we  use  the  low-den¬ 
sity  limit  in  terms  of  Mayer  (  functions 

CjfA, .  A3.  A4 )  .((A, .  Aj  )  /(Aj.  A4  )  /(A4.  A, ) 

(20) 

and  retain  only  the  n  =  1  term  in  the  perturbation  ex¬ 
pansion  18,9) 

/(A,.. 

(1  +.rn(fl,)l  S  \0m]DU.i)r-;\" /r\[  (21) 

II  - 1 

for  the  Mayer  (  function,  where /|-|  (r,, )  is  the  corre¬ 
sponding  Mayer  lunction  for  the  reference  hard-sphere 
system.  The  bridge  diagram  is  now  evaluated  analyti¬ 
cally  using  llankcl  transforms  and  we  find  eventually 
that  when  ;  — 

/}<2)*(oo)'-  rliilf  1 )- fivA'  rrp',’ .  (22) 

in  which  the  contribution  from  the  ihiee-particle  cor¬ 
relation  function  appears  as  the  third  virial  coefficient 
for  hard  spheres!  Adding  tins  to  ( 16).  after  using  (9), 
we  have 

*  L  247rp« y  j  sgrrpO 

where  —  emphasizes  that  our  calculation  of  the  bridge 
diagram  of  0(A* )  is  correct  only  to  lowest  order  in 
Pl  and  Wi .  This  is  sufficient  however  to  test  the  mo¬ 
lecular  theory  against  thermodynamics  using  (18)  for 
the  dielectric  constant.  Instead  of  ( 1 7)  we  now  need 

p^8c/dp?~(e  l)"'.3y  f,(6  Dm-,  (24) 

which  is  con.sisteni  with  ( 18)  to  O(y').  W'e  have  thus 
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shown  that  the  bridge  diagrams  make  significant  con¬ 
tributions  to  the  density  when  r  and  it  appears 
likels  that  the>  alfeet  the  density  profiles  also  when 
the  dipr'les  aie  closer  to  the  wall.  Theoretical  studies 
of  the  ordering  of  dipoles  and  ions  at  a  wall  that  have 
appeared  so  lar  ignore  these  bridge  diagrams. 


5.  The  polari/ation  density  /‘(““./'i)  in  an  open 
system 


From  the  delinilKm  ol  /'(r.  )  in  the  grand  ensem¬ 

ble  |b)  we  have 

/»(.-.£,  l^twip',’  '21  /'/),, tr.fs.nilsiinjldfl.  (251 

In  our  discussion  F.  i  is  independent  of  z  [see  (2)1 . 
Using  the  invariant  expansion  for  /is|{z.f2-  >• 

find,  after  doing  the  necessary  angular  integrations, 
that 

f-i  I  =  Wjp'i'a  s|  (.s  cos'd,  -s  1 )'  ■  e,.  (2b) 

which  is  independent  of  any  approximation  for  the 
wall  particle  correlation  function!  Hmploying  the 
OHNC  apprt'Mination  lor  Fv,!  given  in  (12), 

+  i>  (  ‘ 


=  l(f 


X 


1)  dTrll.''/.;,  l-e  +  '<1 

e  l]  ' 
,1(1  -r  Fvf,  I  2c  +  IJ  ■ 


(2b) 


When  A;,  =  0  (i  e,  in  the  absence  of  clectrostriction  as 
in  the  MS  and  LIINC  approximation)  and  when  all  the 
higher  coeflicients  of  l-.’i  (n  >  2)  in  (3)  are  assumed 
zero,  we  recover  the  constitutive  relation  ( 1 ).  Con¬ 
versely  when  K/i  and  the  higher  coefficients  in  (3)  are 
not  zero,  non-linear  terms  in  the  polarization  density 
will  appear.  Substitution  of  ( I S)  and  (9)  in  (28)  follow¬ 
ed  by  expansion  of  the  denominator  draws  out  the 
term  of  0(/:  2  1  in  the  QUNC  approximation  for  the  po¬ 
larization  density  which  is 

^yiiNt  l)/4i7j[3f2/(2e -r  l)| 


X 


(e  l 
2c  el 


+  0(F^. 


(29) 
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If  the  non-linear  terms  in  (3)  are  neglected 
^(.)l INC  I'f  l)'4n|/ 


I  + 


24rp';p- 


(e  1)-' 
2c  +  1 


/- 


(.7(1) 


.-\  more  ci'inplele  calculation  c'f  this  term  in  ihe 
QIINt  appri'ximatiori  would  leipme  ilie  deteimm.i 
tioi.  ol  the  Coefficient  of  the  l  i  term  m  ( .t  I  when  (■  - 
“  0.  .-Xu  exact  calculation  of  the  lerm  ol  0(/  ■'  I  would 
leijuiie  II-  adiliiion  that  the  contribiilions  ot  the  bridge 
eh.igi.ims  to  electiostriction  [see  tZdl]  .iiid  ti'  the  rel.i- 
Tioiishii'  between  A  and  /  ,  be  included. 
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The  adsorption  of  dipoles  at  a  wall  in  the  presence  of  an 
electric  field:  The  RLHNC  approximation^^ 

John  M.  Eggebrecht,  Dennis  J.  Isbister,“  and  Jayendran  C.  Rasaiah 

Oi'purtrnvni  of  Chemistry,  Vniversity  of  Maine,  Orono,  Maine  04469 
(Rctfued  12  May  1980;  accepted  10  June  1980) 

2  hc  adsorption  of  dipolar  hard  spheres  in  the  presence  of  an  external  electric  field  has  previously  been 
studied  within  the  context  of  the  mean  spherical  approximation.  In  order  to  quantify  the  significance  of  the 
physical  trends  found  above,  the  problem  is  solved  within  the  higher  order  closure  rules  afforded  by  the 
linearized  hypcrnetted  chain  approximation.  1  xpressions  for  the  reduced  dipole  moment  and  tlic 
electric  field  strength  are  derived  using  only  the  asymptotic  forms  of  the  direct  correlation  functions  It  is 
found  that  the  favorable  ohentational  correlations  between  the  dipolar  hard  spheres  and  the  wall  are 
underestimated  by  the  mean  spherical  appro.vimalion.  This  is  emphasized  in  the  enhanced  adsorption  of  the 
dipiilar  species  (at  the  wall  ilselO  for  dipoles  oriented  close  to  the  direction  of  the  field.  However,  the 
nonphysical  features  of  the  mean  spherical  approximation  Imanifested  in  the  negativity  of  the  density  profile) 
are  not  fully  rectified  by  the  use  of  the  linearized  hypernetled  chain  approximation 


I.  INTRODUCTION 

The  adsorption  of  dipolar  molecules  at  a  wall,  in  the 
presence  of  an  electric  field,  is  of  intere.st  in  the  study 
of  electrode  and  membrane  phenomena.  Here  the  ad¬ 
sorption  phenomenon  is  delineated  by  the  distribution  of 
molecule.s  at  a  particular  orientation  !2,  and  distance  r 
from  a  hard  planar  wall  p,(;. .  Ej,  f2i).  The  electric  field 
Ej  emanates  from  this  wall,  the  declination  of  the  field 
with  respect  to  the  wall  beinp  allowed  by  its  nonconduc- 
tive  properties.  Isbister  and  Freasier'  have  investi¬ 
gated  thi.s  problem  for  hard  dipolar  spheres  against  a 
hard  wall  using  the  mean  spherical  appro.ximation 
(MSA).  Their  results  for  the  density  profile  (iilf.,  Ej,  fij) 
of  diixiles  are  of  great  interest  even  though  they  suffer 
from  the  defect  that  the  wall  particle  density  profile 
Pi(z.  Ez,  51,)  assumes  negative  value.s  at  certain  relative 
orientations  of  electric  field  (Ej)  and  dipole  moment  m, 
of  the  particles  in  the  fluid  (the  dipole  orientation  is  de¬ 
noted  here  by  fl,).  However,  the  argument  leading  to 
the  electric  field  at  the  wall  is  not  swayed  by  the  appro.x- 
imation  used,  and  may  be  employed  with  more  accurate 
theories  such  as,  for  example,  the  linearized  hypernetted 
chain  (l,HNC)  approximation.  While  these  theories  could 
be  expected  to  produce  better  results,  they  do  suffer  from 
the  necessity  of  employing  numerical  methods  to  a 
greater  extent  than  is  needed  to  determine  the  density 
profiles  in  the  mean  spherical  approximation. 

This  paper  is  devoted  to  a  study  of  the  wall -particle 
density  profile  using  the  linearized  hypernetted  chain 
approximation^  for  the  wall  particle  and  particle-parti¬ 
cle  interactions,  except  that  the  effects  that  are  inde¬ 
pendent  of  the  orientations  of  the  electric  field  and  the 
fluid  dipole  are  treated  exactly.  In  our  study,  these 
are  the  interactions  between  the  hard  cores  in  the 
fluid,  and  also  the  interactions  between  these  cores  and 
the  hard  wall.  We  call  this  the  renormalized  linearized 
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hypernetted  chain  approxim.ition  (RLHNC)  after  the 
nomenclature  introduced  by  Stell  .iiid  Weis.^  Our  re¬ 
sults  for  this  theory  are  .m  improvement  over  the  cor¬ 
responding  M.S  appro.Nimation  when  the  system  is  char¬ 
acterized  by  (a)  a  reduced  fluiddcnsity  p^jj^nOf 
(b)  a  reduced  dipole  moment  n;’  v  '  0.  5 

(or  equivalently  a  reduced  temperature  T*  ^  I  of  2), 
and  (c)  a  reduced  external  iTeciric  field  £*  =  in, 

of  8/3.  In  contrast  to  the  .MS  .ipproximation,  the 
RLHNC  wall-particle  density  lunctions  Pi(z.  Ej,  52,)  arc 
only  marginally  negative  near  the  wall  for  a  dipole  ori¬ 
entation  in  direct  opposition  to  the  field  (see  Fig.  2). 
These  functions,  however,  can  become  negative  over  i 
larger  distance  :  from  the  wall  when  the  reduced  dipok 
moment  is  increased  to  1.0  without  altering  the  reducea 
electric  field  or  reduced  densily  (f  ic.  5).  This  sug¬ 
gests  even  higher  order  terms,  beyond  the  RLHNC  ap¬ 
proximation,  must  be  included  in  Hu  theory  when  the 
dipole  moment  /;/,  and  the  c.Nlornai  electric  field  are 
both  large. 

II.  GENERAL  THEORY 

The  technique  of  producing  a  w.ill  .next  to  a  fluid  bj 
taking  the  limiting  behavior  of  a  binary  mixture  (with 
densities  p,,  pj  and  radii  H,,  Hfl  detailed  by 

Um  lim  (2.1) 

is  well  known.*  Isbister  and  Freasier'  have  extended 
this  to  introduce  an  electric  field,  a.s  well,  by  consid¬ 
ering  the  corresponding  limit  for  a  diixilar  mixture  un¬ 
der  the  restriction  that  the  dipole  moment  of  par¬ 
ticle  2,  which  eventually  becomes  Hie  wall,  divided  by 
the  cube  of  the  radius  of  the  excluded  volume  Rj,  ^Rz  •  R, 
is  a  constant  (hereafter  called  £o)( 

Um  »/j/R|,  =£„  .  (2.2) 

R  g-  * 

In  taking  these  limits,  in  the  specific  order  pj-  0,  Rj 
—  the  volume  of  the  system  is  allowed  to  grow  faster 
than  Rj,  keeping  p,  constant  through  the  constraint  that 
PzRz—  0, 

The  magnitude  and  direction  of  the  electric  field  Ez 
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Ihe  “wall  dipoif."  whu-h  has  ivccdcci  lo  a  distance 
from  the  wall.  At  that  distance,  the  vi-ctor  r^,  hecoine^ 
pcriwndiculai'  to  tlie  wall,  and  the  electric  held  which 
makes  an  ancle  »,  with  respect  te  this  luirm.il  is  rr.  en 
by  the  appnipri.ite  solution  to 


eosa,  Cf  •  ) . 


2  cose 


(3  cos'^t,  .  1) 


(2.  101 


I  K..  1.  I ’u-  eie.  in,,  iK'l.i  ai.Lil,.  ui  tia  watt  shilled  iigains! 
::k  inclinini  in  e  ‘1  hie  ‘Vail  li;  lie"  a'  e-  '  rc'-i-K-.s  m  minus 
lll'Illi:  iilld  .''laxjaees  'lie  electi  ;e  liel-i  i  .  I'la  anyle  S,  is  the 
liiclinalliin  el  a  aisele  in  llR'  is.li.  11,11,1. 


follow,  when  this  hmit  is  .ipphcd  to  the  dipole-dipole 
interaction  eneruv 


The  solution  to  this  euualion  pifes  Die  direction  ol  me 
electric  field  unuiuely.  in  terms  of  the  orient. then  oi 
the  wall  dipole  .it  minus  infinity  (see  hip.  1). 

The  density  profile  Ej.  S2,)  of  dipolar  molecule' 

at  a  dist.inee  c  troiii  tin  w.ill  depends  on  'ihe  dipoli-  'i  i- 
ent.itior.  f!i,  oe  Ihe  liireclion  ol  the  electric  lielu  E;  .me 
on  Ollier  parara  ters  of  Die  system.  It  is  olit.iiiu  ci  li  oi:, 
the  relation' 

E^,  h,)  lim  Inn  ,i,  id,;,  h,,  S.‘,l  -  if.  12.11' 

Kj  -  ‘  I  L'2'  "  -  j 

wliere  ,i,  is  Die  hulk  density  ol  spi  eie.s  1  .ind  Ididp.. 
h,)  is  ihe  toi.i!  eoi'rel.iliun  lunetion  ol  sjni  les  1  .,mi  2 
.1  hin.ii'y  mi.xture.  The  hitter  is  the  .solution  to  me  tirii- 
stein-Zeriiike  rid.  it  inn 


I  1  V  , 

IlSjaCji*  C,.  , 

‘1  ,  =  I 


D(2.  1)  . 

(2.3) 

where  "  !  tlTyilV.^  is  . 
angular  inlegr.ilions. 

llZlil'zi-  f)))  -/'2l(>  Jl 

'*  2  i 

(2.4) 

'■2l(’'21'  •^2-  ^^1^  Zi 

Here 


■s,  .ind  .sj  are  unit  ^-ccUirs  in  the  direction.':  ol  the  dipole 
moment  vectors  ni,  .ind  nij.  respectively,  (Srj,  rj,  -U) 
IS  the  di(xile-dipole  interaction  tensor,  in  which  rjt  is 
a  unit  vector  in  the  direction  .iloai;  the  line  joinini’  par¬ 
ticles  1  and  2,  .ind  U  ns  Die  unit  tensor.  Tlu'  potential 
m.iv  be  writlcn  in  terms  ot  the  electric  field  Ej,  pro¬ 
duced  .It  the  ioe.iiion  ol  p.irtich  1  l>’.  tin-  seeond  parti¬ 
cle  2. 


(2.  12) 


(2.  13) 

?,(r2,)  A(2,  1)  -  cfifi-j, 1/1(2.  n  . 

(2.  HI 

.ind  usini;  Wortheim’s  mulliplic.ilion  table"  in  Fourii'r 
spaee  lor  the  .mpul.ir  inlei.r,ilions.  the  Ornslein-:  i  rnita 
relalion  roduees  to  three  tviu.ilions 


'•'21  "  ''21  '  ^  D,  h;,  • 


(2.  Ifi' 


V  2i(r^j .  f;,.  f;  I )  -  -  11. 1  s  1  •  E;  . 
who  re 

Ez  -  e,  . 
r  t, 

In  Eq.  (2.6' 

02  -  ISi  ji  12,  -  UI  •  Ss-  (Sr,,  eosr-  - 


(2.5) 


(2.6) 


(2.7) 


.  (2.16) 


,1.4* 


(2.  17) 


^^21  ■■  hji  *  o  "y .  ('>(2(1' 

when  •  !  ilTj  IS  a  convolution  involving  only  sp.ni.il  in 

ler.ralions.  In  Eqs.  (2.  16)  .ind  (2.  17) 


and  hj  is  the  ani;le  which  the  diixile  embedded  in  particle 
2  makes  with  i  j,  (.see  Fip.  1).  Since  the  magnitude  of 
e,  IS  (3cos^b2-  1)'  ^ 


('aefr)  c4,,(,) -3  f  Ts.v-'4,(s)  l2.  IH) 

and  its  inverse  is 

3  r^, 


(3  cos^b,  •  !)'■' ’('2  . 


f  - 


12.  19) 


(2.8) 


where  is  a  unit  vector  in  the  direction  ot  Ej.  Changing 
variables  to  i'2i  -  r.  and  taking  the  wall  limit,  the 
electric  field 


with  similar  expressions  lor  c^alr)  and  Cojl)  )  m  which 
o,  d  =  1  or  2.  Inside  the  hard  cores 


*;«(»■) -1  . 


(2.  20) 


Ej  cFoOcos^bj-  l)'''r2  . 


(2.9) 


which  shows  that  E2  is  independent  of  the  distance  z  from 
the  wall,  but  that  its  magnitude  is  determined  by  the 
strength  (through  E^)  and  the  direction  (through  6^  of 


=  ''Sp(»'>  -  0  ,  r-  Rf  . 

so  that  from  Eq.  (2. 18) 

4p(r)  -  -  3/t-,,  .  r  ‘  «,3.  (2.21) 

where 
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Ka,  -  I  h^,(s)s'',ls  .  (2.22) 

•^K^e 

Follow! lit;  We rlhe  1111.  ^  ttie  relations  for /i2i()  )  and 
/ij,(r)  may  be  uncoupled  by  l.ikint;  linear  combinations 

'<;«(>•)  -  [''««(>■'  *  i  1  .  (2. 23) 

Kiir)  -  KeO  )  -  )  I  3A-„,  .  (2.  24) 

with  analogous  expressions  for  e' ^(i )  and  eiefr).  The 
closure  conditions  (2.20)  when  applied  to  Fqs.  (2.23' 
and  (2.24)  are  equivalent  to 

/C(r)--1  .  .  (2.25) 

On  using  the  linear  combiii.itions  in  Fqs.  (2.  16) 

and  (2.  17),  one  finds 

/'Ii(<-)  -  (-;,()•)  .  A„  /,,•  1,1  .  <-U  .  (2.  26) 

where  p *  -  2p,,  p"  — p,.  Fiiu.aions  (2.  15)  and  (2.  26'  re¬ 
semble  the  Crnsteiii'/erinke  rel.iiions  for  a  binary 
fluid,  in  which  the  molecular  interactions  are  spheric.ii- 
ly  symmetrical. 

In  the  limit  Aj-  ,  E'q.  (2.  19)  together  witli  Eqs. 
(2.23)  and  {2.24)  yields 

/ifi(4-)  =  (2/i*2i(.-)-/iii(c)-3lA',,  .  (2.27) 

A-s  pj- 0,  Eqs.  (2.15)  .ind  (2.26)  reduce  to 

/i|,(r)=<1,(i')  -p,;i;,‘  <•;,  (2.28) 


/ij,(r)  -  <'2i(r)  *-  '“n  ■ 

which  can  be  written  in  bipolar  coordinates  as 


mined  In  Die  inler.ictions  between  the  particles  in  the 
bulk  fluid.  An  analogous  equation  can  be  written  (or  the 
angularly  independent  part  of  the  wall-particle  lota)  cor¬ 
relation  function 

K.l.O  -  e^,(c)  .  2iip,|^.  [AV) [/)'(-) -V'{i)\ 

-  P,/r>;|,(  v)lA>)-B*(U- v|)|  I  ,  (2.35) 

where  /!'(.- '  and  A*(r)  have  definitions  which  corresp'ind 
exactly  to  !i‘{:)  and  D‘iz). 

In  Kqs.  (2.34)  and  (2.35).  Ic5,(<').  ct,(s)i  and 
|•‘l^(>),  are  sets  of  two  different,  but  consistent,  direct 
correhition  functions  with  corresponding  closures  Id:  -  - 
cussed  in  tlie  next  section)  for  the  wall-particle  aiui 
particle-particle  interactions,  respectively.  Tlie  furc- 
tio/i.s  |■|^(•’')  .ind  rj,(r)  for  the  bulk  fluid  particles  are  ct)- 
tained  in  .in  independent  calculation  from 

I.'i'ifr)  -<•;,()•)  ^p,/i5i  •  Cn  .  l2.3t.' 

I/ii(r)  -  I'ljfr)  ♦  A'li  p,  Iiij  *  Cj]  .  (2.2( 

which  art  the  one-component  analogs  of  Eqs.  (2.28)  ..no 
12.  29)  first  derived  in  a  seminal  paper  by  Werthf  in;. 
These  equations  carry  their  own  closures  for  //,,(r)()  .  p 

and  (-nl'  Hr  An)-  We  do  not  actually  solve  Eos.  (2.  36' 
.ind  (2.35)  since  nearly  exact  results  are  available  from 
the  work  of  ferlct  and  Weis®  for  hard  spheres  and  from 
the  study  of  Waisman.  Henderson,  and  Lebowitz’  for 
hard  spheres  against  a  wall.  Our  RLHNC  approximation 
(Sec.  Ill'  implies  that  /i|,(r)  is  the  exact  wall-p.irticii. 
tot.il  correlation  function  (or  hard  spheres  against  a 
hard  w.ill.  Equation  (2.37)  has  been  solved  by  Wei - 
tlieiiii.  in  tlie  mean  spherical  approximation.®  whiie 
Paley  and  his  colleagues  have  treated  it  in  the  1  ilN'C 


/i|i()')  -  r|i(r) f  i/I//i2i(/l  /  i/.s  s  r'nfsl  .  .ind  QHN'C  approxiniations, "  using  the  Verlet-Wei 

^  ''■•I'  tl...  I . a  _  U-r,  .....  tl.f. 


(2.30) 

On  substituting  r  -  A.  •  .  I  -  H.  -  v  .ind  on  taking  the  w.dl 
limit  find,  (or  .-  ’0. 

li'ziU)  --  rjifr)  *  2A',,  j  ,/\  I’/Jit  v)  /  ds  (■;,(.*■)  . 

(2.31) 

where  i  is  the  distance  of  the  center  of  the  dipolar  h.ird 
sphere  from  the  wall.  The  integral  between  the  limits 
and  may  be  simplified  by  noting  that  Iij,(  v)  -  -  1. 

when  yi  0,  and  in  addition  - yl  =r  -  v.  when 

c  >0  and  -«  c  y  q.  On  changing  the  order  of  integra¬ 
tion  between  -«>  and  0. 

(  dv/i|i(v)  f  (/.'  >■(•;,(>■'  .  /  d.<i(r  -  .s)  .sc',i(.<)  . 

a—  -'u-yl 

(2.32) 

Defining  the  functions 

B*(^l  =  [‘  r\,{s)  .v  i/.v  ,  /)•(.0  -  (  *  c'n(s)  .  (2.  33) 

ao  -'ll 

the  wall-particle  ±  equations  become 

2itA„p;ic  [A*(»l  -fl*(c)|-  [/;*(‘ic) 

■2-K^^p\  f  r/v //*,,( \)  |A*(^) -S*(;c  -  ly  )|  . 

-'ll 

(2.. 14) 

where  the  (unctions  A’(.  '  and  />'{.-)  are  entirely  deti  r- 


ory  for  tlu-  liard-sphere  interactions.  We  arc  llicia  (eri 
lelt  with  ilie  necessity  of  solving  only  Eq.  (2.34'  undi  r 
tile  .ippropi-i.ite  closures  for  the  wall-particle  .ind  p.ir- 
tich'-p.iriK  le  direction  correlation  funclions. 

The  const.ints  A’j,  and  A'j,  determine  the  dipole  mo- 
nienl  ii.  j  and  the  electric  field  Ej  through  the  relalio';.' 

tr,.'}.'  *  -  (AiZA'iiPiAji)  -  tg,(- A'nPiA'u'  (2.38' 


^  Aji  [2iA(2AiipjAii)-i  t3.(  —  AjiP]  A  j] '  ,  ( 2 .  .19 

where  /,  and  T  arc  the  Boltzmann  constant  and  absolute 
teniperature,  respectively,  and  V,.(pj  Aj,)  is  defined 
by; 

<,'..iOiA,,)  l-4rpi  I  ef,(r,  p,)  r^i/r  .  (2.40) 

-'o 

Kqu.ition  (2.38)  has  been  derived  by  Wertheiii  in  the 
me.in  spherieal  approximation,  but  its  extension  to  t:.' 

1  IINC  lor  HEHNC)  approximation  is  straighllorw.iid; 
nevertheless  we  present  it  for  completeness  and  as  .i 
prelude  to  the  derivation  of  (2.39),  whleh  Ishister  .uui 
1  re.isier'  disvussed  in  the  mean  spherieal  .ipproxiii..: 
Inin.  The  derivation  of  Eq.  (2.38)  rcst.s  on  llu  .mi  r:,;i 
u>li(  fiii  M'  of  '•u(i  )  |.s<’e  Kqs.  (3.1)  and  (3.6' 
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rf,() '  -  l)\2.  1)  ,  .1^  I  -  ■  (2.41) 

On  insrrtini;  tins  in 

,'i’i(r)  -cfiir)  -2  I  ,/>  '-■‘.'[■(o,  ,  (2.42) 

one  sees  that  is  siioi’i  aiui  lends  to  zero  as 

r  —  : 

asr-'  .  (2.43) 

When  llie  inverse  relation  ^ei.  Fcis.  (2.  19)  and  (2.18)| 


Alue!:  ult-nlieu!  to  Kij.  (2.39).  11  .-should  he  eiiiph.j- 

-a/ed  that  I  q  .  {2.3b)  aiui  (2.39)  are  relalions  dei  i\- 
Irojii  tin  a>>in|)loiie  lorin.s  ol  ej,  and  e-,  in  the  Ms  .la  ; 
I.IINC’  .ipprnMnialnnr- .  l  ln*  dijxde  nionu-nt  //. ;  .iml  i;.- 
eleeirie  lieUl  tv  .iri*  determined  Iroin  A'^  .md  Aj,,  a!o 
the  sulutinns  to  tq.  (2.37)  havt-  been  obtained. 

The  enii-t.iiu  A'21  also  retrieves  llie  eoi  ll u  ient.s 
ami  ),  .Much  appear  in  the  invariant  expansion  t,l 

the  wall -parlieli-  correlation  liiiuTions 

E;.  //;,(.  )  •  //?,(.-)  yt(2.  1)  -  a.{2.  1)  (2. 


<-n('-erf,('  '  -  4  (  r-., (si  (2.44) 

)  -V 

i.s  ta)ten  to  tiu'  hnut  i  -  '  ,  oiu'  .il.sn  tuids  tli.U 


ah 

kT 


/' 


(  Jl(s)  S^c/.s 


j  [2ct|(s.  2A',,  ,i.l  ■is,(s.  -  A'liP,)  !n^(/s  . 

(2.451 

wlierc  we  liavo  u.sed  the  an.iiints  i>:  l  iis.  (2.23)  .tnd 
(2.24)  applied  to  rfits)  in  tlie  l.ist  stnp.  Tbe  re.suit 
Stiven  in  Eq.  (2.38)  follows  inimodi.iteU  on  .ipplyinp 
Eq.  (2.40). 


The  relation  between  A’j,  .ind  £  ,  is  .ilso  derived  from 
the  usymptotie  from  ol  cfilr)  in  the  \I8  .tnd  LHNC  ap- 
pro.ximations  [see  Eqs.  (3.1)  and  '3.H)  in  8ec.  in|.  The 
equation  for  dipolar  mixtures  wliieli  eorresponds  to  Eq. 
(2.45)  is® 


I'lm 
'21  •• 


111^)1/1  3A'2i  (■ 

,.5*,,r  =  -  lull  I 

IjlAi  '21'“  »  21  •'(> 


[2<-2,(s)  '  rjil.slj.arfs  , 


(2.46) 

where  we  have  not  canceled  the  1*21  in  the  denominator, 
because  we  intend  to  take  the  wall  limit.  Substituting 
t'21-^2'-'.  and  taking  the  limit ,  reduces  the  left 
hand  side  of  Eq.  (2.46)  to  kT .  .ind  replaces  the 

upper  limit  of  the  integral  in  Eq.  (2.46)  by  .  The 
Fourier  transform  of  Eq.  (2.29),  wliicli  i.s  an  Ornstcin- 
Zernike  equation  for  mixture.s,  .i.s  .u  ~  0,  yields 


?5.(A0=[1-A'„p;  r!,(£)|/;^,(/o) 


(2.47) 


from  which,  when  £  =  0,  we  have 


r  Caif.s)  ,s®rfs  = 
•'o 


[i-x,.pi  ?;,(o)|  I 


/l|,(,s)  . 


(2.48) 

In  the  above,  the  “'"represents  the  three  dimensional 
Fourier  transform 


/(/?)  =  y^rfrexpOk- r)/(iri ) 


and  fsc  or  h. 

Since  /ij,(.s)  is  a  short-ranged  function  equal  to  -  1  for 
s<Rt,,  it  readily  follows  that 

lim  T  f  fili(.s)  -  5  .  (2. 49) 

R  2*  “  r 2j  Jq 

Hence. 


s  Ajj  (2  [  1  -  /f  11  Pi  rii(O)  I  •  1 1  -  A 11  Pi  Cii(O)  I  I 


(2.50) 


by  iiiverlinc  !  <|.  (2.23)  and  (2.24)  .ind  using  Eq.  (2.27-. 
Both  ;£.(.  )  (for  tlie  closure  rules  considered  in  Die  m 
section)  .Old  TfA:)  .ire  shurt-r.inged  functions  .ind  It  na 
to  /ero  .IS  .  .iiid  the  .isymplotic  (orm  of  i' 

thereloi  i 

lim  i'll'B  I  3A ji  .  (2.  72 

.cliicl'.  loiiows  ir'nii  Eq.  (2.27)  .md  the  Diet  Dial  I'D' 

.liso  short  r.i'i.,cd.  Tlie  .isymplotie  lorm  ol  Du  w.,11- 
p.iriieir  coin  l.itioii  function  lor  the  BLHN'C  cio.-uri  . 
disiussed  in  Du  next  section,  coincides  with  th.i!  -  ;  D,. 
M.S  api)ro\i:i..itioM'  .iiid  is  Diei  elore  given  by 


Bin  E;,  3A-,,ii(2,  1) 

(2.  53' 

3ii,  ,  £,/i(2.  1) 

K’  I  _  2(^  2A  ]  t  f  A 

'2.  54' 

3/';i^(2.  i)| 

'  r*  i2v^2K;,,).K^,)  • 

(2.5'' 

where 

(2.  5t 

and 

J'*  '  ^1' 

//.  1  (3  cos^bg  “  1)^'  ‘  n/j 

(2.  .57' 

Since  Die  electric  field  is  independent  of  the  distance 
from  the  wall.  £2.  fl,)  does  not  decrease  to  zero 

as  ,  unless  the  magnitude  of  the  electric  field  01 

V{2,  1)  is  also  zero.  The  angular  average  of  /i2i(.',  E  , 
is  however  zero  in  the  limit  z-"-': 

lim  /  /i2i(.'.  Ej,  n,)rff2|  =  0 

*-•  J 

(2.56> 

since  the  angular  averages  of  D{2,  1)  and  .i(2,  1) 
and 

are  zero 

lim  7i|,(z)  =  0  . 

(2.  59) 

We  shall  now  consider  the  details  of  these  equations  m 
the  context  of  the  mean  spherical  and  linearized  hyper - 
netted  chain  closure  rules. 

HI.  THE  CLOSURE  RELATIONS 

The  closures  for  the  wall-particle  direct  correl.ition 
functions  arc  readily  derived  by  taking  the  wall  limit  oi 
the  ctirro.s|xmding  closures  in  the  bulk  fluid. 


J.  Chem.  Phys.,  Vol.  73,  No.  8.  tSOctotipr  1980 


.!98-l 


Eggebrecht,  Isbister,  and  Rasaiah:  Dipoles  at  a  wall 


A.  MSA 

Till'  I'losurt'  tor  tlu'  ilirei't  rri'l.itioii  function  is 


Tin-  .itioic  i'<|U.itiiins  (3.121  'Acre  solved  itii  .dr.s  i  . 
1  oiiij.  :ii  r  Ui  ooiicr.itc  Ilic  ..ccoiiiji.irn  inn  ficuri  , 


lor;-/;,,,  (3.1) 

kl  >■ 

ivliicli  i.s  i'i)uiv.ilcnl  to 

(■;,^()  i  0  ,  for  )■  ■  .  (3-  2) 

This  i.s  un.iflictiHl.  ior  IT,,  on  l.ikiim  the  wall  lunit 
ot  iT|( r'i . 


B.  LHNC  approximation 

startinc  with  the  h\  pi  rnctti  d  chain  (UNO  appro.\inia- 


r„,(2,  1)  1)  -  lni,h,..(2,  n  -  f,a(2,  1)  *•/■  (r  •«„,) 

(3.3) 

ami  usin'^  the  invariant  expansions  of  1)  and 

/)„a(2.1)  :cf.  Eqs.  (2.131  and  (2.  14)  |.  one  find.s 


kr 


/'(2.  1)  . 


where  is  the  spherically  symmetric  part  of  the  po¬ 
tential.  E.xpandinp  the  luparithm  up  to  the  linear  term, 
.ind  coUectiii!;  and  comparinv  coefficients  of  1.  Di2,  1), 
and  .i(2, 1).  we  have 

■-  I'iiir)  -  InA'ialr)  -  V>  a(i  )  kT  ,  (3.  5) 


.■Sa(v)0,Sa(r)[l-A'l.(.-)-'i‘^  • 


IV.  RESULTS  AND  DISCUSSION 

Tin  siuatmn  of  Eq.  (2.51).  lor  .  E  ..  1. ;  ■  ■  ’ 

t.iiiied  .1.-  tile  confluence  of  three  distinct  coiijjr.ii.i tiom 


(II  The  hulk  correlation  functions  ivi  .  c.ilcui.do  c. 
lilt  in.iiuier  dt  scribed  by  Patey  ,  from  Eg.  (2.37  .  J  , 
reduced  dilx.le  moment  dependence  ol  this  rel.de.r;  , 
(.lined  in  till- constant  fv  1,1  throunh  Eg.  (2.33:. 
lie  e.xprissed  in  llie  MSA  and  ULHNC  appru.xin..dio:i- 


wliere  and  (>n(/fi*i)  are  the  value,-  ol  r;';irl  .me. 

h(.(i  )  ininiediatcly  outside  contact.  The  value  o!  ,T.; 
■v.i-  .idjusled  until  this  difference  assumed  the  desiri’' 
reduced  riiixife  moment.  The  iteralion  of  Eg.  (2.37 
w.is  perlornied  usinp  t.isl  Eourier  transforir  lecr.nni;..  - 
.iiid  mixed  solutions  to  spei-d  convcrtienee.  Tin  ouh, 
tluid  correl.ition  funrlion.s  obtained  showed  i.wi 
.e_reetiie.'it  iiith  those  of  Patey." 


(2'  The  electric  field  which  eini  rites  Iron,  tin 
w. IS  determined  by  Eqs.  (2.9),  (2.39).  and  (2.40!.  J;.: 

field  ancle  a,  of  Fit;.  1  was  taken,  in  separate  i  .ili  .il,.- 
tions.  .IS  O'.  45',  and  90".  The  field  angle  is  iel..U‘c. 
to  tile  wall-dipole  orientation  t/j  through  Eq.  (2.  lOi.  I'le 
solution  of  which  appears  in  Fig.  3.  The  magnitude  ul 
the  reduced  electric  field 


E,;/??,  ^(3  cos«bj-lV^^ 


c^s(v)=ltS2(r)[l  -,gi«(r)-V  . 

(3.7) 

Defining 

(3.8) 

bU(r)=/,UrHi-jrlBirr''  . 

(3.9) 

with  b^eSr)  and  ()„,(»-)  also  defined  by  equations  analogous 
to  Eqs.  (2.18),  (2.23),  and  (2.24),  respectively,  the 
last  two  equations  can  be  written  in  the  form 

cls(r)  =6*j(r)  ,  .  (3.10) 

In  the  RLHNC,  Eq.  (3.5)  is  replaced  by  the  exact  closure 
for  hard  spheres  against  a  wall. 

In  the  wall  limit,  df,(e)  =*21(2),  and  making  use  of 
Eq.  (2.27),  the  closure  condition  for  the  wall -particle 
direct  correlation  function  becomes 

c‘^,(^)  =  [lx/^|l(^)][l-l?|.(^)''),  ^>0  .  (3.11) 

On  substituting  this  in  Eq.  (2.34),  we  have,  for  2>0, 

h\,{z)  =  [Aii(z)  -  1 K  ifli(«)  2itK„  p\  [B*(»)  -  B‘(*)  I 

v[B‘(«)-^*(2)K  /  f4(£)lB*(*)-B*(U-,vl)Uv|  • 

(3. 12) 

When  ,g|,(<-:)  =1,  the  formal  equation  for  h\^{z)  in  the 
MSA  approximation  is  recovered. 


w.is  t.iken,  as  in  the  earlier  work  of  Isbister  an, 
Frcasier,  to  be  with  the  reduced  dipole  momcn* 
squared  m''  fixed  at  either  0.5  or  1.0.  For  purposi. 
of  comp.irison  of  these  parameters  to  a  molecular  sy-- 
um,  the  reduced  di|x>le  moment  squared  for  HCl  (e  ; 

.  1.03  D)  .It  275 'K.  assuming  a  diameter  of  3.5  A.  i- 
approximately  0.65.  A  reduced  field  strength  01  8  3 
for  this  system  corresponds  to  an  electric  field  of  near¬ 
ly  1.9  •  10^  V '  m  or  a  surface  charge  density  a  of  1  elec¬ 
tronic  charge '1000  A^. 

Wfc  have  also  carried  out  calculations  at  the  same 
surface  charge  density  (or  electric  field  E^)  for  a  fluid 
at  a  reduced  density  of  0.  7  with  the  reduced  dipole  mo¬ 
ment  III*  -  2.0.  These  numbers  correspond  approxi¬ 
mately  to  those  appropriate  for  liquid  water  (m,  s  1. 85 
D,  /?,,  =  2.  76  A)  at  room  temperature.  The  reduced 
electric  field  [which  contains  m,  and  bJ,  in  its  defini¬ 
tion  (4.2))  is  now  only  0.71.  (It  may  be  useful  for  the 
reader  to  bear  in  mind  that  it  is  an  artifact  of  our  defini¬ 
tion  of  El  that  an  increase  in  the  dipole  moment  m,  re¬ 
sults  either  in  a  reduction  of  BJ  when  the  surface  charge 
density  is  held  constant,  or  an  increase  in  the  surface 
charge  density  when  eJ  is  unchanged. ) 

(3)  The  spherically  symmetric  part  of  Eq.  (2.51). 
i.  c. ,  /;2i(.'),  was  determined  using  the  technique  of 
Waisnian,  Henderson,  and  Lebowitz.  ’  The  function 
/rjifr)  wa.s  first  obtained  from  Eq.  (2.28)  as  the  Percus- 
Yevick  solution  and  then  corrected  to  produce  an  essen- 
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6  -1 


^  t- 


K,  5R, 


KK'».  Ilie  waU-pariiole  distriljution  functions  ^  j  *2,  E  .  ii-\ 
as  a  function  of  the  distaive  j  from  the  uall  for  tlifferenl  oritn' 
tallons  iJ.  of  the  fluid  dipoles  wliere  £7  =  s  3,  p.* 

-  u.oTls  and  a.  -  ii  .  A.  is  the  radius  of  a  fluid  dipole.  - - 

lU.H Nt'  apprnximui ion.  -----  MS  approximation. 


-2- 


7 


b  t', 


5P, 


I  I'J,  Wall-partu'lf  ills' riiiution  funcUons  for  the  system 
(lepieieil  in  fig.  exeepi  tiu.i  o,  i-.j.j'. 


FIG.  ■).  Wall-particle  distriliullon  functions  for  the  systems 
depicted  In  Figs.  2  and  2  except  that  a,  =  90°. 


3R 


FIG.  ■>.  ITie  wall-particle  distribution  tunciion  .c  .ic,  E.,  ii,i 
for  the  system  depicted  in  Fig.  I  except  that  in'-  l.o. 


120-  ■■  “  ^ 

I 


0  PC  P 


I  K..  The  expansion  coefficients  *2)  and  ict  :ts  n  lune- 
lion  01  the  wall-particle  distance  a  when  £?  =  H  :i  and  ,i?  -  o,  ;i7;i. 
the  upper  and  lower  curves,  in  each  case,  are  lor  o  f  i.p 
ando..i,  respectively. 
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IK'i.  7.  I'he  \v;  U-p:ir:ioK' luiu'tioii  ir 
as  a  luiiL'lion  ol  liu'  tlislanss  It  ■'  !;,<.■  wall.  U»r  ditlori-nl 

orU'iilalions  1.’,  o!  Itif  iUai  -  when  E*  ‘  ."I,  tn*  l.u, 

pT  U.V,  tin^i  t';  I'iii’  la-  i'iva-il  paiMnioU-i's  nt*  ;tn<l  .»*  arf 

tbo.so  appi'opriaU'  lor  litiui'i  w.tior  a  room  u-mperaturo 
-l.S.'t  1),  '^1]  -  7. 76  liu-  I'Oihaavl  (.•loctric  liolU  ft  .0.71 

corresponds  to  a  surlaco  densitv  td  1  ckadronic  charite  liiOO  \  . 


tiaUy  exact  wall-pariiclo  correlalion  (unction  in  com¬ 
parison  to  Monte  Carlo  simulation.  The  reduced  densil.v 
of  hard  .spheres,  upon  which  depends  solely,  was 

taken  as  p*  -PiKii  =0.  573. 

The  computations  were  comhined  through  Eqs.  (2.23). 
(2.  241,  and  (3. 12)  to  yield  the  coefficients  of  the  angular 
functions  of  Kq.  (2.51).  The  integrals  over  bulk  corre¬ 
lation  functions  (2.33)  as  well  .is  the  integral  of  Eq. 

(3.  12)  were  performed  by  trape/.oidal  approximation. 

The  Iterative  solution  to  Eq.  (3.12)  was  rapidly  conver¬ 
gent.  Finally,  the  expansion  eoelficients  were  inserted 
into  Eq.  (2.  51). 

In  Figs.  2-5  we  present  our  results  for  .i,'j,(/,  Ej,  fi,) 
in  the  KLHNC  approximation.  The  expansion  coeffi¬ 
cients  h?,(r)  and  /(2i(r)  are  sliown  in  Fig.  6.  A  compari¬ 
son  of  RLHNC  and  MSA  results,  for  a  field  angle  of  zero 
and  reduced  dipole  moment  squared  of  0.5  and  1.0,  are 
shown  in  Figs.  2  and  5,  respectively. 

In  the  case  of  w*^  =  0.5.  E*  -  8  3.  and  p*  -0.573 
fiiiU,  Ej.  U,)  in  the  RLHNC  approximation  assumes 
negative  values  only  very  near  and  at  contact  for  a  di¬ 
pole  orientation  in  opposition  to  the  field.  For  higher 
values  of  wt^  =  1.0,  at  the  same  reduced  electric  field 
and  fluid  density,  these  regions  are  more  pronounced 
and  extended  in  range  (Fig.  5).  An  increase  in 
from  0.5  to  1.0  at  constant  E*  and  pj,  however,  corre¬ 
sponds  to  a  fourfold  increase  in  the  surface  charge  den¬ 
sity  front  1  electronic  charge  1000  to  one  every 


250  Holli  the  RLHNC  and  M.SA  treatments  result 
in  distribulions  showing  exhaneed  adsorption  lor  f.ivor- 
able  diiKile  orientations,  but  tiu'  repulsive  inter. o  tion- 
of  the  eli’etrie  field  with  unfavorably  aligned  dijioles  i- 
more  iTe.irly  visible  in  the  RLHNC  approxiiualion.  I'h.^ 
theory,  like  the  MSA.  is  essentially  linear  in  cli.ir.ietei  . 
and  cannot  prevent  the  distribution  functions  from  bf  - 
eoniing  iieg.iiive  when  the  dipoles  are  aligned  .ig.tinst  li.i 
ficUl.  if  .ii  the  same  time  the  theory  predicts  .i  l.irgi  er- 
haneemeiu  ol  ibe  .idsorplion  of  dipoles  aligned  with  llie 
field.  When,  tor  example,  the  eleetrie  field  is  perja  i;- 
dieul.ir  to  llie  wall.  botJi  theories  ))rediet  that  the  den- 
sil\  proliles  isee  Figs.  2  .ind  5)  are  syminetricai  .diou' 
llie  prolile  U.r  diixiles  perpendieul.ir  to  the  held  f:  . 

-  t:  2  or  3:  21.  A  large  cnh.ineeiiu-Ml  of  dijiole.s  .ilp:.')' . 
with  (he  Held  would  lead  to  an  equally  large  depleimi;  , 
dipoles  orienled  against  the  field,  which  m.i\  req.nri 
negative  w.ill-p.irticle  dislributitm  lunctions.  ?vir;hi  : 
improvements  beyond  the  RLHNC  .ipproximalion  so  ih: 
then  be  neeessarv.  Figure  7  shows,  however.  Ui.it  ll.’ 
Hl.HNC  theory  provides  plausible  density  prolio  s  .  v.  :, 
when  the  squ.ire  ol  the  dipole  moment  n.  *'  m  inei  v.is,  ,! 
to  4.0  hut  the  surface  charge  density  i.s  maint.iinei.  .it  1 
eleetjonie  charge' 1000  A“. 
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'll.  '.  IsbisI,'!'  and  li.  i‘.  I'reasler.  .1.  Stai.  I’fivs.  2e.  ! 

Jhe  noialion  lhai  i.s  ii.sed  in  this  paper  Is  shgl-.li' 
leren;.  In  partionlar,  /?„  is  the  radius  ol  the  n  spe.'U> 
ratlu-r  than  the  tlianieler.  and  £.  and  E,  .as  deliiied  ;ier,  .art 
h  ami  £  r.  respect ivel.v,  of  Ishister  .ami  I  reasier. 

The  l.ilXt  approximation  lor  diiiolar  lluuls  propostxl  1.'  ii.  '  . 
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